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ABSTRACT 

’The  present  study  concerns  the  use  of  an  analytic  model  for  closed 
brain  injury  to  evaluate  General  Motors'  Severity  Index  for  head  injury. 
This  index,  which  considers  the  relative  importance  of  the  duration  and 
intensity  of  the  pulse  by  means  of  an  impulse- integration  procedure,  has 
been  demonstrated  empirically  to  be  an  excellent  means  of  comparing  the 
severity  of  pulses  of  varying  shapes.  Four  pulses  (square,  half-sine, 
skewed  and  triangular),  each  having  the  same  Severity  Index  but  different 
magnitudes,  were  employed  as  inputs  to  the  model. 

The  model  response  to  each  of  the  four  impacts  was  determined. 
Graphs  of  the  time  variation  of  stresses  in  the  shell  and  pressures  in 
the  fluid  at  both  the  impact  pole  and  at  the  counter  pole  demonstrate 
excellent  correlation  between  model  response  and  Severity  Index, 
independent  of  the  shape  and  magnitude  of  the  impact  force. 


Criteria  to  assess  the  severity  of  impacts  classically  have  been 
based  upon  specific  parameters  of  the  pulse  wave-form.  Maximum 
force  or  acceleration  sustained,  rate  of  onset,  pulse  shape,  duration, 
and  the  rate  of  change  of  acceleration  are  representative  of  character¬ 
istics  of  the  pulse  which  have  been  selected  as  injury  criteria.  However, 
as  von  Gierke  *  has  noted,  the  situations  in  which  any  single  aspect  of 
the  impact  wave-form  characterises  the  system  response  are  limited. 
This  observation  is  substantiated  by  the  work  of  Kornhauser  and  Gold^ 
who  have  found  that  the  incidence  of  injury  depends  on  both  the  magnitude 
and  the  duration  of  the  impact.  Their  tolerance  curves  of  velocity 
change  (Atr)  versus  average  acceleration  (a)  indicate  that  threshold 
values  for  both  Av  and  a  must  be  exceeded  before  injury  will  occur. 

Thus,  a  relation  between  both  magnitude  and  duration  must  be  considered 
in  assessing  impact  severity. 

Based  on  the  assumption  that  the  severity  of  an  impulse  wave¬ 
form  is  a  function  of  both  duration  and  intensity,  Gadd^  developed  an 
exponentially  weighted- impulse  criterion  for  estimating  head  injury 
resulting  from  frontal  impacts.  From  the  work  of  J .issner4’  5  in  which 
tolerance  curves  of  intracranial  pressure  versus  time  indicated  that  the 
time  required  to  produce  severe  concussion  in  dogs  decreased  as  the 
pressure  increased,  and  from  that  of  Eiband^  which  demonstrated  a 
similar  trend  for  sled-test  tolerance  curves  of  acceleration  versus 
time,  Gadd  noted  that  both  results  demonstrated  a  downward  sloping 
tolerance  curve  in  the  range  of  vehicle  occupant  cranial  impacts.  He 
further  noted  that  when  the  data  was  represented  on  a  log-log  plot  the 
curve  could  be  approximated  by  a  straight  line  at  an  angle  of  less  than 
45°  with  the  horizontal.  He  concluded  that  an  impulse- area  criterion 
would  not  satisfactorily  predict  injury.  The  resulting  severity  index 
which  Gadd  devised  is  depicted  in  Figure  1. 

In  general,  the  severity  index  (S.  I. )  may  be  expressed  mathemati¬ 
cally  as: 


5.1.  -j (a)n  dt 


where  S.  I.  is  the  severity  index  (a  number);  a  is  the  acceleration, 
pressure,  or  force  producing  an  injury  threshold  of  a  specified  degree; 
n  is  the  exponential  weighting  factor;  and  t  is  time  in  seconds. 


For  head  injury  resulting  from  frontal  impacts,  a  is  the  accelera¬ 
tion  of  the  head  expressed  in  "g"  units,  and  n,  based  upon  empirical 
results,  is  2.  5.  The  exponential  weighting  factor  signifies  the  impor¬ 
tance  of  the  high  intensity  segments  and  the  insignificance  of  the  low 
intensity  portions  of  the  impacts,  while  the  integral  accounts  for  the 
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effect  of  the  duration  of  both  the  high  intensity  and  the  low  intensity 
portions  of  the  (rurve. 

l 

The  purpose  of  this  study  is  to  investigate  the  usefulness  of  the 
severity  index  in  comparing  the  severity  of  pulses  of  different  peak 
magnitude  and  varying  wave-form.  The  severity  index  has  been  used 
primarily  to  investigate  the  injury  potential  of  various  pulses  of  one  to 
fifty  millisecond^  duration  applied  to  the  frontal  bone  of  the  head.  A 
previously  developed  mathematical  mo^el  for  head  injury^  has  been  used 
therefore  to  study  various  pulses  having  the  same  severity  index  but 
differing  in  wave-jform  and  peak  magnitude. 

The  head  injury  model  assumes  tliat  the  skull  is  a  linear,  homo¬ 
geneous,  isotropi^,  thin  spherical  shell ,and  that  the  brain  is  an  ideal 
(acoustic)  fluid,  losing  extensional  shell  theory  and  assuming  axisymmetry 
of  the  load,  the  following  set  of  three  coupled,  second  order,  linear 
partial  differential!  equations  with  variable  coefficients  obtain  in 
spherical  polar  coordinates:  '< 
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Figure  2  depicts  the  relationship  between  the  displacements  u,  v,  and 
w  and  the  coordinate  axes. 

Equations  (1)  -  (3)  are  solved  by  finite  difference  methods  using 
backwards  difference  analogs  and  employing  techniques  described  in 
references  7  and  8.  Figure  3  shows  the  segmentation  of  the  shell  for 
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Figure  2.  Coordinate  System  used  to  Describe  the  Model. 


Figure  3.  Shell  segmentation  for  Finite  Difference  Analogs. 
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application  of  the  finite  difference  analogs. 

The  model  may  be  applied  to  an  input  with  arbitrary  temporal  and 
spatial  distributions,  provided  axisymmetry  is  maintained.  Four 
different  inputs  were  selected  to  evaluate  the  severity  index,  and  the 
model  responses  for  each  have  been  compared.  The  four  acceleration 
pulse  shapes  chosen  to  represent  those  which  might  be  encountered  in 
actual  vehicle-interior  impact  situations  are  square  pulse,  half- sine 
pulse,  triangular  pulse  and  an  arbitrary  impact  pulse  as  suggested  by 
Goldsmith^.  These  four  acceleration  pulses  are  depicted  schematically 
in  Figure  4.  The  mathematical  representations  of  the  four  pulses  are 
shown  in  Figure  5,  in  which  t  is  the  independent  variable,  time,  and  T 
is  the  duration  of  the  impact. 

A  pulse  duration  of  0.  0001  sec  was  selected  for  each  of  the  four 
input  functions.  The  intensity  of  the  pulses  was  selected  by  requiring 
that  each  of  the  four  pulses  yields  the  same  severity  index.  That  is, 

(a).2,  5  dt'  =  500. 

The  input  to  the  model,  represented  by  the  term  Fext  in  equation 
(2),  has  units  of  lbf/in2.  Each  of  the  four  input  pulses  was  expressed 
in  these  units  by  assuming  a  head  weight  of  10  pounds  and  an  even 
distribution  of  the  load  over  a  polar  cap  of  15°  polar  angle.  Figure  6 
depicts  the  four  inputs  as  functions  of  both  time  and  polar  angle. 
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Numerical  Results 


The  following  properties  of  the  shell  and  the  encased  fluid  are 
used  in  the  calculations: 


Shell: 


_  ,  ,n6  lbf 2 

E  =  2x10  /in 


H  =  0,  25 
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Ps 


=  3  in. 

=  0.  15  in. 
=  2.00x10 


lb  sec2..  4 
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Figure  4.  Schematic  Representation  of  the  Four  Input  Pulses. 
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Figure  5.  Mathematical  Description  of  the  Four  Input  Pulses. 
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Fluid: 


K  =  300,  000  lbf/in2 
p£  =  0. 938x10  /m 

C  =  56,  553  in/sec. 
o 


Using  this  data  and  determining  the  model  response  to  each  of 
the  four  input  pulses  yields  the  results  shown  in  Figures  7-11.  The 
resting  position  is  considered  to  be  a  zero  datum,  and  all  pressures 
are  gage  rather  than  absolute  pressures. 

Pressure  variation  with  time  at  polar  angles  of  3°  and  177°,  the 
positions  approximating  the  impact  pole  and  the  counter  pole,  respec¬ 
tively,  are  shown  in  Figures  7  and  8.  Fluid  pressures  are  studied  at 
r  =  1.  0  (immediately  adjacent  to  the  shell)  since  this  location  is  where 
the  greatest  pressure  fluctuations  occur.  The  fluid  at  the  impact  pole 
is  compressed  initially  and  as  the  impact  load  decreases  in  intensity  the 
compression  diminishes  and  the  fluid  is  subjected  to  a  tensile  stress. 

The  duration  of  this  tensile  stress  is  approximately  125  ^sec.  The 
pressure  at  the  counterpole  does  not  reflect  a  change  from  the  zero 
datum  until  the  stress  wave  in  the  shell,  which  generates  instantaneous 
centers  for  wave  propagation  as  it  traverses  the  shell,  approaches  the 
counterpole.  As  the  shell  attempts  to  separate  from  the  fluid,  tensile 
stresses  occur  in  the  fluid.  These  stresses  are  more  diffuse,  less 
severe  and  longer  lasting  than  the  tensile  stresses  at  the  impact  pole. 
Although  the  solution  has  not  been  obtained  for  a  period  that  is  sufficiently 
long  to  determine  the  duration  accurately,  it  is  estimated  that  the  tensile 
stress  at  the  counterpole  is  maintained  for  approximately  300  ijsec. 

The  variation  of  the  shell  radial  displacement  as  a  function  of 
time  for  both  the  impact  pole  and  the  counterpole  are  shown  in  Figures 
9  and  10.  The  load  produces  an  initial  indentation  at  the  impact  pole 
followed  by  an  elastic  recovery  as  the  load  is  removed.  There  is  a 
continual  outward  displacement  at  the  counterpole  which  assists  in 
producing  the  tensile  stresses  in  the  fluid  at  that  point.  Membrane 
stresses  in  the  shell  are  shown  in  Figure  11  for  the  impact  pole. 

Initially  high  compressive  stresses  are  generated  during  the  impact 
period,  followed  by  elastic  recovery  when  the  load  is  relieved. 

Each  of  Figures  7-11  shows  the  model  response  for  the  four 
impact  pulses  selected.  Inspection  of  these  five  graphs  reveals  a 
similar  model  response  for  the  four  pulses  although  the  pulse  shapes, 
rates  of  onset,  and  peak  intensities  differ. 
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NONDIMENSIONAL  TIME,  t 

Figure  7.  Pressure  versus  time,  0  =  3°,  r  s  1.0 


Figure  8. 


NONDIMENSIONAL  TIME,  t 

/SZ- 

Pressure  versus  time,  0  =  177  ,  r  =  1.0 
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NONDIMENSIONAL  TIME,  t 

Flyure  11.  Normal  stress  in  the  6-direction  versus  time,  <b  -  3° 

/V 

135 


i/ 


Conclusions 


The  intent  of  this  study  has  not  been  to  determine  the  effectiveness 
of  the  General  Motors'  Severity  Index  in  predicting  injury  thresholds, 
since  this  result  could  be  accomplished  only  in  a  program  which  combined 
analytic  and  experimental  or  clinical  investigation.  The  objective  of 
this  study  has  been  to  assess  the  value  of  the  Severity  Index  in  predicting 
a  single  parameter  equivalent  for  pulse  wave-forms  of  varying  shapes 
and  magnitudes.  As  evidenced  by  the  results  shown  in  Figures  7-11, 
the  Severity  Index  does  appear  to  be  an  excellent  means  of  comparing 
impact  functions. 

The  current  work  may  be  extended  to  an  investigation  of  model 
response  to  pulses  of  durations  within  the  range  of  known  vehicle  interior 
impacts,  or  to  a  study  of  pulses  of  identical  peak  intensities  but 
different  durations.  The  predictive  capability  of  the  severity  index 
under  these  conditions  could  therefore  be  determined. 


Acknowledgments 

The  author  expresses  his  appreciation  to  the  Life  Sciences 
Division  of  Technology  Incorporated  for  support  of  this  study. 


136 


NOMENCLATURE 


a 

a 


h 

K 


n 

P  (0*  t) 
P' 

P 

R 

SI 

r' 

r 

T 

t' 

t 


Radius  of  the  shell  middle  surface  (L) 

Average  acceleration 

Constants,  determine  magnitude  of  the  acceleration 
functions 

Constant,  determines  decrement  rate  of  the  loading  function 

Compressional  wave  speed  in  the  encased  fluid,  f  K/pQJ*^(L/T) 

Wave  speed  in  the  shell,  [E/p^l  -  m2)]1/2(L/T) 

Young's  modulus  for  the  shell  material  (F/L2) 

.  2 

External  loading  function  (F/L  ) 

Shell  thickness  (L) 

Bulk  modulus  of  the  fluid  (F/L2) 

Exponential  weighting  factor 

2 

General  loading  function  (F/L  ) 

Pressure  in  the  fluid  measured  with  the  undisturbed  state  as 
aero  datum  (F/L^) 

Nondimens ional  pressure,  P'/K 

Value  of  the  discrete  radial  variable  at  the  shell  boundary 
Severity  index 

Dimensional  radial  coordinate  (L) 

Nondimensional  radial  coordinate,  r'/a 
Nondimensional  duration  of  loading 
Dimensional  time  (T) 

Nondimensional  time,  C  t'/a 


u' 

u 

V* 

V 

w' 

w 


x,y,e 


Av 

e 

T1 


0,0 


\ 

H 

Pf 

Pi 

♦ 

<t> 

Subscripts 

h 

i 


Dimensional  meridional  displacement  (L) 
Nondimens ional  meridional  displacement,  u'/a 
Dimensional  parallel  displacement  (L) 
Nondimensional  parallel  displacement,  v'/a 
Dimensional  radial  displacement  (L) 
Nondimensional  radial  displacement,  w'  /a 
Cartesian  coordinates 
Incremental  change  in  velocity 
Dimensionless  parameter,  a/h 
Dimensionless  parameter,  ap^/hp^ 

Spherical  coordinates 

Dimensionless  parameter,  Co/Cs,  the  speed  ratio 

Poisson's  ratio 

Density  of  the  fluid  (FT2/L*) 

Density  of  the  shell  (FT2/L*) 

Dimensional  velocity  potential  (L2/T) 
Nondimensional  velocity  potential,  f/aCo 

Refers  to  discrete  variables  in  the  radial  direction 
Refers  to  discrete  variables  in  the  0  direction 
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